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Abstract 

The no-shp boundary condition at a soHd-Hquid interface is at the center of our understanding of fluid me- 
chanics. However, this condition is an assumption that cannot be derived from first principles and could, in 
theory, be violated. In this chapter, we present a review of recent experimental, numerical and theoretical 
investigations on the subject. The physical picture that emerges is that of a complex behavior at a liquid/solid 
interface, involving an interplay of many physico-chemical parameters, including wetting, shear rate, pressure, 
surface charge, surface roughness, impurities and dissolved gas. 
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1 Introduction 



The vast majority of problems in the dynamics of Newtonian fluids are concerned with solving, in particular 
settings, the Navier-Stokes equations for incompressible flow 

p(at + u • V)u = -Vp + /iV^u, V-u = 0. (1) 

The list of problems for which this task has proven to be difficult is long. However, most of these studies 
assume the validity of the no-slip boundary condition, i.e., that all three components of the fluid velocity on a 
solid surface are equal to the respective velocity components of the surface. It is only recently that controlled 
experiments, generally with typical dimensions microns or smaller, have demonstrated an apparent violation 
of the no-slip boundary condition for the flow of Newtonian liquids near a solid surface. 

We present in this chapter a tentative summary of what is known about the breakdown of the no-slip 
condition for Newtonian liquids and discuss methods and results of experiments, simulations and theoretical 
models. This topic is of fundamental physical interest and has potential practical consequences in many areas 
of engineering and applied sciences where liquids interact with small-scale systems |147l I149| , including flow 
in porous media, microfluidics, friction studies and biological fluids. Furthermore, since viscous flows are 
relevant to the study of other physical phenomenon, such as the hydrophobic attraction in water, a change 
in the boundary condition would have significant quantitative impact on the interpretation of experimental 
results [T71[T751[T77| . 

The present chapter complements previous work [631 11531 1176) and is organized as follows. In ^ we 
present a brief history of the no-slip boundary condition for Newtonian fluids, introduce some terminology, 
and discuss cases where the phenomenon of slip (more appropriately, this may often be "apparent slip" ) has 
been observed. In ^ we present the different experimental methods that have been used to probe slip in 
Newtonian liquids and summarize their results in the form of tables. A short presentation of the principle and 
results of Molecular Dynamics simulations is provided in [21 weW as remarks about the relation between 
simulations and experiments. We then present in fJJlan interpretation of experimental and simulation results 
in light of both molecular and continuum models, organized according to the parameters upon which slip has 
been found to depend. We conclude in offering a brief perspective on the subject. 

2 History of the no-slip condition 

2.1 The previous centuries 

The nature of boundary conditions in hydrodynamics was widely debated in the 19th century and the reader 
is referred to |^ for historical reviews. Many of the great names in fluid dynamics have expressed an 
opinion on the subject at some point during their careers, including D. Bernoulli, Euler, Coulomb, Darcy, 
Navier, Helmoltz, Poisson, Poiseuille, Stokes, Hagen, Couette, Maxwell, Prandtl and Taylor. In his 1823 
treatise on the movement of fluids jll5| . Navier introduced the linear boundary condition (also proposed later 
by Maxwell |ll)9p . which remains the standard characterization of slip used today: the component of the fluid 
velocity tangent to the surface, , is proportional to the rate of strain, (or shear rate) at the surface,^ 

U|| = An- (Vu+ (Vu)^) • (1 - nn), (2) 

where n denotes the normal to the surface, directed into the liquid. The velocity component normal to the 
surface is naturally zero as mass cannot penetrate an impermeable solid surface, u • n = 0. In Eq. A 
has the unit of a length, and is referred to as the slip length. For a pure shear flow, A can be interpreted as 
the flctitious distance below the surface where the no-slip boundary condition would be satisfled (see Fig. 
Note that on a curved surface the rate of strain tensor is different from the normal derivative of the tangential 
component of the flow so all the terms in Eq. (0) need to be considered ^S! • 

A century of agreement between experimental results in liquids and theories derived assuming the no-slip 
boundary condition (i.e., A = 0) had the consequence that today many textbooks of fluid dynamics fail to 

^Alternatively, the right-hand side of this boundary condition, when multiplied by the shear viscosity of the liquid, states that 
the tangential component of the surface velocity is proportional to the surface shear stress. 
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A = < X < oo X = oo 

Figure 1: Intepretation of the (Maxwell-Navier) slip length A. 

mention that the no-slip boundary condition remains an assumption. A few monographs however discuss the 
topic. In his classic book (p. 576), Lamb realizes that no-slip is the most probable answer but leaves the 
possibility open for extraordinary cases: 

It appears probable that in all ordinary cases there is no motion, relative to the solid, of the fluid imme- 
diately in contact with it. The contrary supposition would imply an infinitely greater resistance to the 
sliding of one portion of the fluid past another than to the sliding of the fluid over a solid. 

Similarly, Batchelor |S] (p. 149) offers two paragraphs where the question is discussed in detail, including the 
role of molecular effects in smoothing out discontinuities. He also mentions the importance of an experimental 
validation of the no-slip condition: 

The validity of the no-slip boundary condition at a fluid-solid interface was debated for some years during 
the last century, there being some doubt about whether molecular interactions at such an interface lead 
to momentum transfer of the same nature as that at a surface in the interior of a fluid; but the absence 
of slip at a rigid wall is now amply conflrmed by direct observations and by the correctness of its many 
consequences under normal conditions. 



2.2 Terminology 

We introduce here some useful terminology that is used throughout this chapter. 

• Phenomenon of slip: Refers to any situation in the dynamics of fluids where the value of the tangential 
component of the velocity appears to be different from that of the solid surface immediately in contact 
with it. 

• Molecular slip (also intrinsic slip): Refers to the possibility of using hydrodynamics to force liquid 
molecules to slip against solid molecules. Such a concept necessarily involves large forces |15H| . Let 
us denote by ct a typical molecular length scale and by A the Hamaker constant for the intermolecular 
forces. Molecular slip will occur when intermolecular interactions 0{A/a) are balanced by viscous forces 
0{fj,a'^j) where fi is the shear viscosity of the liquid and 7 the shear rate; this can only happen for a 
very large shear rates 7 « A/fxa^ ~ 10^^ s~^, where we have taken the viscosity of water /i = 10^'^ Pa-s, 
and typical values A « 10~^^ J and a « 0.3 nm. 

• Apparent slip: Refers to the case where there is a separation between a small length scale a where the 
no-slip condition is valid and a large length scale L ^ a where the no-slip condition appears to not be 
valid. Well-known examples of such apparent slip include electrokinetics [I'SS^ (in this case a is the the 
thickness of the double layer) and acoustic streaming ;8, (in this case a is the thickness of the oscillatory 
boundary layer). Similarly, a liquid flowing over a gas layer displays apparent slip (see ^jSJ. 

• Eflfective slip: Refers to the case where molecular or apparent slip is estimated by averaging an appro- 
priate measurement over the length scale of an experimental apparatus. 
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2.3 Traditional situations where slip occurs 

The phenomenon of shp has already been encountered in three different contexts. 



Gas flow. Gas flow in devices with dimensions that are on the order of the mean free path of the gas 
molecules shows significant slip |113| . An estimate of the mean free path is given by the ideal gas formula 
~ 1/ (v^Ti'O'^p) where p is the gas density (here taken as the number of molecules per unit volume); for air 
under standard conditions of temperature and pressure, im ~ 100 nm and, in general, im depends strongly 
on pressure and temperature. The possibility of gas slip was first introduced by Maxwell He considered 

the flow of an ideal gas and assumed that a percentage (1 — p) of the wall collisions were specular whereas a 
percentage (p) were diffusive. Such an assumption allows an exchange of momentum between the gas and the 
wall. The corresponding slip length is given by 

^ ^ 2(2 -p) 

The case of a rough surface with only specular reflections was considered in ^3] . In general, a Knudsen number 
defined as the ratio of the mean free path to the system size Kn = im/ L is used to characterize the boundary 
condition for gas flow, with slip being important when Kn ^0.1 j5()l . 



Non-Newtonian fluids. The flows of non-Newtonian fluids such as polymer solutions show signiflcant 
apparent slip in a variety of situations, some of which can lead to slip-induced instabilities. This is a topic 
with a long history and is of tremendous practical importance, and we refer to Refs. |2(J[ I42L W7\ 1751 1^ ll(J2l 
II 1(JI 11421 and references therein for an appropriate treatment. 



Contact line motion. In the context of Newtonian liquids, molecular slip has been used as a way to remove 
singularities arising in the motion of contact lines, as reviewed in |43[ I169j . Solving the equations of motion 
with a no-slip boundary condition in the neighborhood of a moving contact line leads to the conclusion that 
the viscous stresses and the rate of energy dissipation have non-integrable singularities. It was first suggested 
[75| that a local slip boundary condition for the flow would remove the singularity, and indeed it does [7D1 1168| 
(see also |17l)p . Furthermore, since the slip length appears via a logarithmic factor in a condition involving 
the contact line speed, it has virtually no influence on macroscopic quantities such as force and pressure drops 
on scales larger than the capillary length = {l/pgY^'^ |158| . where 7 is the liquid surface tension, p the 
liquid density and g gravity. Consequently, macroscopic measurements on moving droplets cannot be used 
in general to deduce the exact slip law gSj. The slip length could also become velocity dependent due to a 
combination of microscopic roughness and contact angle hysteresis '83'. Such local slip near contact lines was 
confirmed by early Molecular Dynamics simulations .94i . In a different context but with a similar purpose, 
slip was used to remove the singularity in the mobility of particles sliding near solid surfaces [411 



2.4 Newtonian liquids: No-slip? Slip? 

The development of the Surface Force Apparatus in the 1970's I 79II81|[T^ (see H3.1|l has allowed for more than 
thirty years of precise probing down to the nanometer scale of both structure and dynamics of many Newtonian 
liquids against mica EHl E3 EHl IZ21 EZl EDI El EOl EH HHHSSl ■ Experimental methods have included squeeze 
and/or shear flow for a variety of polar and non-polar liquids displaying a wide range of wetting conditions 
and shear rates. With the exception of the flow of toluene over Cgo (Fullerene)-coated mica 25 , these studies 
have conflrmed the validity of the no-slip boundary condition and the bulk rheological behavior down to a 
few nanometers. At smaller length scales, an increase in viscous resistance has been reported, with qualitative 
differences between the behavior of water and other non-polar liquids 57, 62, 89, 13(|. The conclusions have 
been confirmed by Molecular Dynamics simulations |135| and are consistent with studies of flow in capillaries 
with diameters of tens of nanometers [3 1^ . 

In this context, the large number of recent published experiments reporting some form of (apparent) slip 
with A ~ 1 nm-1 /xm in the flow of Newtonian liquids is surprising P P HHI [T gl H gj PUl |P |^ |^ |^ | ^ |^ 
EZIEHllMlIinainiinSlinailSniinaiinniinilina aUowed to re-discover a few early 
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Schnell ITITI 
Churaev 



Kisolova ^ 
Cheng ESI 



Cheikh |^ 
Choi |32l 



Surfaces 


Liquids 


Wetting 


Roughness 


Shear rates 


Slip length 


L/NL 


Glass+DDS 


Water 


_ 


_ 


10^ - 10'' s^i 


1-10 fj.ni 


L 


Quartz+TMS 


Water 


70 - 90° 




1 s-i 


30 nm 


NL 




Mercury 


115 - 130° 




10^ - lO'' s-i 


70 nm 


NL 




CCI4 


Complete 


_ 




no-slip 






Benzene 


Complete 


— 


- 


no-slip 


- 


Quartz+CTA(+) 


CTAB solutions 


70° 




10^ - 10^ s^i 


10 nm 


L 


Glass+photorcsist 


Water 




5 A (pp) 


10= - 10" s-i 


no-slip 






Hexane 








10 nm 


L 




Hexadecanc 








25 nm 


L 




Decane 








15 nm 


L 




Silicon Oil 








20 nm 


L 


Poly(carbonatc)+PVP 


SDS solutions 


< 90° 




- lO'^ s"^ 


20 nm 


L 


Silicon 


Water 


a 0° 


11 A (rms) 


10^ - 10^ s"^ 


0-10 nm 


NL 


Silicon+OTS 


Water 


> 90° 


3 A (rms) 




5 — 35 nm 


NL 



Table 1: Summary of slip results for pressure drop versus flow rate experiments. The following symbols 
are used in this table: — : unknown parameter; DDS: dimethyldichlorosilane; TMS: trimethylchlorosilane; 
CTAB / CTA(-I-): cetyltrimethyl ammonium bromide; PVP: polyvinylpyridine; OTS: octadecyltrichlorosilane; 
CCI4: tetrachloromethane; SDS: sodium dodecyl sulfate; pp: peak to peak; rms: root mean square; L: slip 
independent of shear rate; NL: shear rate dependent. 



studies reporting some degree of slip |21l 1341 1^ 11411 1160| . In part, this chapter is an attempt to describe and 
interpret these more recent experimental results. 



3 Experimental methods 

As will be discussed below, a large variability exists in the results of slip experiments so it is important to 
first consider the different experimental methods used to measure slip, directly or indirectly. In these setups, 
surface conditions may usually be modified by polymer or surfactant adsorption or by chemical modification. 
Two broad classes of experimental approaches have been used so far, indirect methods and local methods. 



3.1 Indirect methods 

Indirect methods assume Eq. |(2J) to hold everywhere in a particular configuration and infer A by measuring 
a macroscopic quantity. Such methods report therefore effective slip lengths, and they have been the most 
popular so far. If the effective slip length is A, then a system size L at least comparable L ^ X is necessary in 
order for slip to have a measurable impact. 



Pressure drop versus flow rate. This standard technique is used in many studies [M)\ 1321 IH^ 1881 1141j . 

where the main results are summarized in Table A known pressure drop Ap is applied between the two 
ends of a capillary or a microchannel and the flow rate Q = J u dS is measured. A slip boundary condition 
leads to a flow rate, Q{X), larger then the no-slip one, Qns, by a factor that varies with the ratio of the slip 
length to the system size; e.g., for a circular pipe of radius a, we get 

m^.^^A (4) 

Using this method, we also note that two groups have reported a larger resistance than that expected with 
the no-slip condition in microchannels |122j and for flow through small orifices jHS]- Their results are not well 
understood but might be due to electrokinetic effects or flow instabilities. 

^The dependence of these results on the size of the system was studied in |101| . where there was some evidence that A increased 
with the size of the system. 
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Surf&ccs 


Licjuids 


Wetting 


Rough lie ss 


Slicar rates 


Slip length 




o: rSocniiKe |iDJ 


Silica 


Propanediol 


ft; u 




i s 


no-slip 








PropajnGdiol"4^ Va 








1 /im 






„ 


PDMS 








no-siip 






Silica+DETMDS 


Prop&iicdiol 


70 — 80° 






no-slip 








Propa.nediol-|- Va 














„ 


PDMS 








no-slip 




FR: Pit ll:^5HT26ll68l 


S&phire 


Hcxadccanc 


Complete 


4 A (rms) 


102 _ lo* s-1 


175 nm 


L 




Saphirc+FDS 




65° 






no-slip 






Saphire+OTS 




40° 






400 nm 


L 




Saphire+STA 




25° 






350 nm 


L 


PIV: Trctheway [TS^ITCT 


Glass 


Water 


* 0° 




10^ 


no-slip 






Glass+OTS 




120° 


2 A 




0.9 fim 




PIV: Joseph 113 


Glass 


Water 


a 0° 


5 A(rins) 


10^ s-i 


50 nm 




Glass+OTS 




95° 






no-slip 






Glass+CDOS 




95° 






50 nm 




SP: Churaev [^H 


Quartz 


KCl solutions 




2 nm (pp) 


10= s-i 


no-slip 






Quartz+TMS 


KCl solutions 


80 - 90° 


25 nm (pp) 




5 — 8 nm 




FC: Lumma [TTOl 


Mica 


Water 




15 nm (pp) 


10 s""^ 


0.5 - 0.86 fim 




Glass 


Water 


5 - 10° 






0.6-1 fini 








NaCl solutions 








0.2 - 0.6 fim 





Table 2: Summary of alternative experimental methods to infer slip. The symbols used in this table are given in 
TableQ] with additional symbols as: S: sedimentation; FR: fluorescence recovery; PIV: particle image velocime- 
try; SP: streaming potential; FC: fluorescence cross-correlations; DETMDS: diethyltetramethyldisilazan; FDS: 
perfluorodecanetrichlorosilane; STA: stearic acid (octadecanoic acid); CDOS: chlorodimethyloctylsilane; Va: 
vacuum; PDMS: polydimethylsiloxane; KCl: potassium chloride; NaCl: sodium chloride. 

Drainage versus viscous force. This technique consists in imposing the motion (steady or oscillatory) of 
a curved body perpendicular to a solid surface, and measuring the instantaneous resistive force, which may be 
compared with that from a model of the fluid motion in the gap, assuming no-slip or slip boundary conditions 
|64[I121II131| . This method is similar in principle to the pressure drop vs. flow rate method, with the difference 
that, here, the pressure and velocity fields are unsteady. The two most common narrow-gap geometries are 
either a sphere of radius a close to a planar surface or two crossed cylinders of radius a. For both cases, the 
viscous force F opposing the motion has the form |172j 

F^^r'-^, (5) 

where V is the instantaneous velocity of the moving body, D the minimum distance between the two surfaces, 
and /* the slip factor. If the no-slip boundary condition is valid, f* — l, otherwise when there is slip, /* < 1 
and depends on the slip lengths on both surfaces. The calculation for /* = /sUp in the case of equal slip 
lengths is given by (B^ 

f ^ 

J slip 



D\f 6A , 



(6) 



3A 

and has been extended to account for two different slip lengths |172| and for the case of any curved bodies 
|173j . Note that when D X, /sUp goes to zero as /sUp D\n{6X/D)/3X, so that the viscous force, Eq. jSJ), 
only depends logarithmically on D] this is a well-known result in the lubrication limit j59| . 

Two different experimental apparati have been used to measure drainage forces, the Surface Force Appa- 
ratus (SFA) and the Atomic Force Microscope (AFM). The SFA was invented to measure non-retarded van 
der Waals forces through a gas, with either a static or dynamic method |?mil54| . and was extended in [221 to 
measure forces between solid surfaces submerged in liquids. More recently it has been used by many groups 
to measure slip in liquids, with results summarized in Table |31 This technique usually uses interferometry 
to report the separation distance between the smooth surfaces. The moving surface is attached to a spring 
system of known properties so the difference between imposed and observed motions allows a calculation of 
the instantaneous force acting on the surfaces. 

The AFM was invented by Binnig, Quate and Gerber and has also been used for many investigations 
of slip, with experimental results summarized in Table 01 A flexible cantilever beam (typically, microns wide 
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Surfaces 


Liquids 


Wetting 


Roughness 


Shear rates 


Slip length 


L/NL 


Chan 


Mica 


OMCTS 






10 - 10'^ s"^ 


no-slip 








Tetradecane 


_ 


_ 




no-slip 


_ 






Hexadecane 


_ 


_ 




no-slip 


_ 


Israelachvili 


Mica 


Water 






10 — 10" s^"' 


no-slip 






Tetradecane 


_ 


_ 




no-slip 


_ 


Horn 


Silica 


NaCl solutions 


45° 


5 A(av) 


10 — 10^ s^^ 


no-slip 




(j60rg6S ^^^J 


6 surfaces ^^^^ 


9 liquids 




0.2 - 50 nm (pp) 


1 — 10 s^^ 


no-slip 




Baudry |^ 


Cobalt 


Glycerol 


20 — 60° 




1 - 10"' s""- 


no-slip 






Gold+thiol 




90° 






40 nm 


L 


V^OLLlIl- JjlZOIlIlC l^^l 


'■jiass 


Glycerol 





Inm (pp) 


1 10" s~ ' 


no-slip 






t^iass-|~vj ± 


Glycerol 


0^° 

ifO 






ou — zuu nm 


T 

Li 






Water 








DU — iiuu nm 




Zhu 11921 


AT' 1 TTT^ A 

Mica+HDA 


Tetradecane 


12° 


J5i 1 A (rms) 


10 — 10^ s 


0—1 f-Lm 


NL 


Mica +U 1 cj 


Tetradecane 


44° 






— 1.5 f-LUi 


NL 






Watei 


1 1 no 






n o . . 
U — Z.O fJiUl 


NL 


7 Vin ^^^^ 

zjnu 


iViica-|-w ± 


vvatei 


75 — 105° 


6 nm (rms) 


5 -1 
10 — 10 s 


no-slip 








Tetradecane 


35° 


6 nm (rms) 




no-slip 






Mica+.8 PPO 


Water 


85 — 110° 


3.5 nm (rms) 


„ 




NL 






Tetradecane 


21 — 38° 


3.5 nm (rms) 




— 5 nm 


NL 




Mica+.2 PPO 


W^ater 


90 — 110° 


2 nm (rms) 


„ 


— 20 nm 


NL 






Tetradecane 




2 nm (rms) 




0-20 nm 


NL 




Mica+OTE 


Water 


110° 


0.2 nm (rms) 




- 40 nm 


NL 






Tetradecane 


38° 


0.2 nm (rms) 




- 40 nm 


NL 


Zhu [T^ 


Mica+PVP/PB 


Tetradecane 




ai 1 nm (th) 


10 - 10° s^^ 


no-slip 






Mica+PVA 


Water 






10 - 10° s"^ 


0-80 nm 


NL 


Zhu 11^51 


Mica 


n- Alkanes 


Complete 




10 ~ 10° s^"' 


no-slip 






Mica+HDA 


Octane 






10 - 10° s^i 


- 2 nm 


NL 






Dodadecane 








- 10 nm 


NL 






Tetradecane 


12° 






0-15 nm 


NL 


Cottin-Bizonne l38l 


Glass 


Dodecane 


a ° 


1 nm (pp) 


10^ - 10" s-i 


no-slip 








Water 


a 0° 






no-slip 






Glass+OTS 


Dodecane 








no-slip 








Water 


105° 






20 nm 


L 



Table 3: Summary of slip results for experiments using the Surface Force Apparatus (SFA). The symbols 
used in this table are given in Tables ^ and [21 with additional symbols as: HDA: 1-hexadecylamine; OTE: 
octadecyltriethoxysilane; PPO: polysytrene (PS) and polyvinylpyridine (PVP), followed by coating of OTE; 
PVP/PB: polyvinylpyridine and polybutadiene; PVA: polyvinylalcohol; OMCTS: octamethylcyclotetrasilox- 
ane; av: average; th: polymer thickness. Note that many entries in this table, including the largest slip lengths, 
are from the same group (S. Granick, U. Illinois). 

and hundreds of microns long) with a small (tens of microns) attached colloidal sphere is driven close to a 
surface, either at its resonance frequency or at fixed velocity, and the deflections of the beam are measured. 
Since the mechanical properties of the beam are known, deflection measurements can be used to infer the 
instantaneous drainage force on the colloidal particle. 

Sedimentation. This experimental method was used in 16^, with their results summarized in TableEl The 
sedimentation speed under gravity of spherical particles of radius a is measured. If the particles are small 
enough, their motion will occur at small Reynolds number; in that case, the sedimentation velocity with a slip 
length A, u(A), is larger than its no-slip counterpart, unS: according to 

^ 1 + 3A/a 

t;NS l + 2A/a' ^ ' 

Streaming potential. This is the experimental technique employed in [33], with their results summarized 
in Table[2| A pressure drop is applied to an electrolyte solution between the two ends of a capillary and creates 
a net flow. Since the surfaces of the capillary acquire in general a net charge in contact with the electrolyte, 
the net pressure-driven flow creates an advection-of-charges current which results in a surplus of ions on one 
end of the capillary, and a deficit in the other end. If the two ends of the capillary are not short-circuited, 
a net steady-state potential difference, termed the streaming potential, exists between the two ends of the 



7 



Craig 0Ol 
Bonaccurso 1181 

Sun ITsnl 

Bonaccurso 1171 

Neto ITTfil 

Vinogradova 11781 

Henry ^3 

Cho EU 



Surfaces 


Liquids 


Wetting 


Roughness 


Shear rates 


Slip length 


L/NL 


Silica+gold+tliiols 


Sucrose sol. 


40 - 70° 


6 A (rms) 


10 - 10*^ s"^ 


0-15 nm 


NL 


Mica/glass 


NaCl solutions 


Complete 


1 nm (rms) 


10^ - 10*^ s"^ 


8-9 nm 


L 


Mica/glass 


1-propanol 


< 90° 


Inm (rms) 


10^ - 10*^ s"^ 


10 - 14 nm 




Silicon / Glass 


Sucrose sol. 


Complete 


7 A (rms) 


10^ — 10^ s^^ 


— 40 nm 


NL 


Silicon/Glass+KOH 






4 nm (rms) 




80 nm 


NL 








12,1 nm (rms) 




100 — 175 nm 


NL 


o 1 ii c a ~p go 1 d "p "C 11 1 o is 


oucrose soi. 


4Q _ 7Q° 


fi k { rm«~l 

U r\. 1 i llio J 


10 — 10^ s~ "'^ 


n 1 S Ti m 


NL 


Silica/ glass 


NaCl solutions 


Complete 


o AyTins ) 


10 — 10^ s "'" 


no-slip 




Polystyrene 


!N^aC> 1 solutions 


90° 


2 . 5 nm (rms) 
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Table 4: Summary of slip results for experiments using the Atomic Force Microscope (AFM colloidal probe). 
The symbols used in this table are given in Tables HT^ with additional symbols as: KOH: potassium hydroxide; 
HTS: hexadecyltrichlorosilane. 



capillary and is such that the current due to advection of net charge near the solid surfaces is balanced by the 
conduction counter-current in the bulk of the electrolyte 1741 11321 . If the fluid experiences slip at the 
wall (and if the ^-potential is unchanged by the treatment of the surface), a larger current will occur, hence a 
large potential difference Ay(A) given by 

§W.1 + A.. (8, 

IWnS 

where k is the Debye screening parameter, which gives the typical distance close to the surface where there is 

1/2 _____ 

a net charge density in the hquid, = (ereofcsT'/2e^no) |138| . Here Cr is the dielectric constant of the 
liquid, eo the permitivity of vacuum, ks Boltzmann's constant, T the temperature, e the electron charge and 
no the number density of ions in the bulk of the solution. 

3.2 Local methods 

All of the methods above have the disadvantage that the slip boundary condition, Eq. (O, was not veri- 
fied directly, but instead was estimated via the assumed effect of slip on some other measured macroscopic 
parameters. A few techniques have been introduced that try to alleviate this indirect estimation of slip. 

Particle image velocimetry (PIV). This method was proposed to investigate slip in |87[ 11621 1163| : we 
have summarized their results in Table |21 (see also jSHl 11881 Q,89 ) . Let us consider a pressure-driven flow 
between two parallel plates with separation distance 2h. In this case, a non-zero slip length A leads to a 
velocity field 



TT I \ ^ 



7? 2\' 



(9) 



which shows that a change in the condition at the boundary has a bulk effect: the (no-slip) Poiseuille flow is 
augmented by a plug flow. The idea of PIV is then to use small particles as passive tracers in the flow, to 
measure the velocities of the particles with an optical method and check whether the velocities extrapolate to 
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Table 5: Summary of MD simulation results for Lennard- Jones liquids with N liquid atoms. List of symbols: 
HML: heavy-mass lattice; NN: fixed atoms of carbon nanotube; FL: fixed lattice; BF: flow driven by a body 
force; CF: Couette flow; CL: contact line. 

zero at the solid surface. Since small particles have large diffusivities, results need to be averaged to extract 
the advective part of the tracer motion. Also, the particles in general move relative to the fluid owing to 
hydrodynamic interactions 164', so care is needed when interpreting the measured velocities. 

Near-field laser velocimetry using fluorescence recovery. This is the experimental technique proposed 
in jfiSI 11251 \12(\\ . and we have summarized their results in Table |21 In this method, the velocity field of small 
fluorescent probes is measured close to a nearby surface. An intense laser illuminates the probes and renders 
them non fluorescent (photobleaching). Monitoring the fluorescence intensity in time using evanescent optical 
waves allows to obtain an estimate of the slip length. Note that the fluorescence intensity evolves in time due 
to both convection (part which depends on slip) and molecular diffusion, so a careful analysis is needed. Also, 
because of the fast diffusion of molecular probes, the method is effectively averaging over a diffusion length 
1 /im, which is much larger than the evanescent wavelength. 

Fluorescence cross-correlations. This is the latest experimental method, proposed by 1U8' , and their 
results are also summarized in Table|5J Fluorescent probes excited by two similar laser foci are monitored in two 
small sample volumes separated by a short distance. Cross-correlation of the fluorescence intensity fluctuations 
due to probes entering and leaving the observation windows allows to determine both the flow direction and 
intensity. The measured velocities are averaged over the focal size of microscope and the characteristics of the 
excitation laser. 

4 Molecular dynamics simulations 

Molecular Dynamics (MD) simulations are useful theoretical tools in the study of liquids |21 1^21 1135| which 
have been extensively used to probe boundary conditions. We summarize below the principle of the technique 
and discuss the interpretation of their results in the continuum limit. 

4.1 Principle 

MD simulations integrate numerically Newton's law of motion for single atoms (or molecules) 

j 
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where rrii is the atomic mass, the position of atom i, and the interatomic (or intermolecular) force 



between atoms i and j, that is Fj 



-S/iVij where Vij is the interaction potential. Potentials used in 



simulations range from the Lennard-Jones two-body potential 



V,, 



(11) 



where e is an energy scale, <t the atomic size, and r.y the distance between atoms i and j, to more realistic 
potentials including many-body or orientation-dependent interactions pi I92L [T^ . The set of Eqs. H1U|) are 
integrated in time, with appropriate numerical cut-offs, and with specified boundary conditions and initial 
conditions. Usually initial positions are random and initial velocities are taken from a Boltzmann distribution. 
It is also possible to modify Eq. H10(l slightly to model evolution at constant temperature either by coupling 
the system of atoms to a heat bath or by a proper rescaling of the velocities at each time step. Interactions 
with a solid can occur by adding different wall atoms, either fixed on a lattice or coupled to a lattice with a 
large spring constant, to allow momentum transfer from the liquid but prevent melting. The constants (c^) 
in Eq. I|ll|) allow variation of the relative intermolecular attraction between liquids and solids, which therefore 
mimicks wetting behavior. Using a simple additive model j78|. the case of a partially wetting fluid with contact 
angle 9c can be modeled with 

cos0c = -l + 2^^ (12) 
Plcll 

where ps (pl) is the solid (liquid) density and cls a-nd cll are, respectively, the liquid-solid and liquid- liquid 
intermolecular constants. Finally, two types of flow can be driven. In the first kind, atoms that constitute the 
wall(s) are driven at a constant velocity and the bulk liquid has a Couette-flow profile. In the second kind, 
each liquid atom is acted upon by a body force and the liquid has a Poiseuille-flow profile. 



4.2 Results 

The method described above has been used to study slip in different types of liquids jJl I35L 1361 1371 152L I55L 1661 
IH2»i24i 95, 114, 145, 151 , 156., .157. .158) . with results summarized in Table[Sl Early simulations showed no-slip 
except near contact lines |^ I156| . More recent investigations have reported that molecular slip increases with 
decreasing liquid-solid interactions |7l l35[[Tni . liquid density |95II157| . density of the wall 1158', and decreases 
with pressure ,7 . The model for the solid wall, the wall-fluid commensurability and the molecular roughness 
were also found to strongly influence slip [3 I55L I82L I151L I158| . Related investigations include contact line 
motion [S3 (and references therein) , motion of a sliding plate [0^1 and the validity of Stokes drag formulae at 
small scales |171| . 



4.3 Interpretation in the continuum limit 

Results from MD simulations can sometimes be difficult to interpret in the continuum limit. First, for com- 
putational reasons, simulations to date are limited to tens of thousands of atoms, which restricts the size 
and time scale of the simulated physical system. The three control parameters in the simulations are the 
molecular /atomic mass to, the interaction energy e, and the molecular /atomic size a. Consequently, lengths 
are measured in units of tr (~ 3 A) and times in units of the molecular time scale r '-^ yjma'^ /e (~ 10~^^ s). 
Simulated systems are therefore limited to tens of nanometers, and times scales to nanoseconds. The conse- 
quence of this observation is that MD simulations always probe systems with much higher shear rates than 
any experimental setup. For example, in MD simulations of Couette flow, the typical wall velocity is C/ ~ ct/t, 
corresponding to typical shear rates 7 ~ cr/r/i where h is the typical length scale of the simulation box, usually 
a few tens of a. Consequently, 7 ~ 10^^ s~^, which is orders of magnitude larger than experimental shear rates. 
Note that this does not apply to investigations inferring slip length from equilibrium simulations |14[ I15| . 

A second signiflcant issue in interpreting results of MD simulation was pointed out by Brenner and Ganesan 
in the case of particle diffusion near a solid surface |19| . and concerns the scale separation between molecular 
and continuum phenomena. The idea is that the correct boundary condition in the continuum realm should 
arise asymptotically as a matching procedure between the outer limit of the inner (molecular) system and the 
inner limit of the outer (continuum) system. By doing so, the change in the physical behavior within a few 



10 



intermolecular length scales of the surface is explicitly taken into account, which allows to make a distinction 
between conditions at a boundary and boundary conditions. As a consequence, slip lengths should not be 
measured literally at the molecular scale but arise as the extrapolation, at the boundaries, of the far field 
hydrodynamic results, a procedure which is not always performed appropriately. 

5 Discussion: Dependence on physical parameters 

Having described the different methods by which slip is investigated, we present in this section a discussion 
of both experimental and simulation results and compare them with theoretical models. The discussion is 
organized according to the physical parameters upon which slip has been found to depend. 

5.1 Surface roughness 

Roughness influences resistance. Be it at the molecular size "KF or on larger scales [171 1581 1?^ 11251 1194j . 

roughness and geometrical features have been observed to influence the behavior at liquid-solid interfaces. Not 
only does roughness leads to an ambiguity as to the exact location of the surface, but it impacts the dynamics 
of the nearby fluid, leading experimentally either to an increase [581 182l[l25lll94| or a decrease ^3 of friction 
with roughness. 

Roughness decreases slip. The physical idea for a roughness-induced resistance is straightforward: on 
the roughness length scale, a flow is induced that dissipates mechanical energy and therefore resists motion. 
For the same reason, a bubble with a local no-shear surface rises at a finite velocity in a liquid. More 
generally, geometrical features of size a on a surface can be solely responsible for a large resistance on large 
scales L ^ a, independently of the details of the local boundary condition on the scale a. This feature was 
first recognized by Richardson |134| who assumed a periodic perfectly-slipping surface shape and performed 
asymptotic calculations for the limit a/ L — > 0; in this limit the no-slip boundary condition was recovered (see 
also |117[I118] ). The calculation was revisited by Jansons 8^ who considered a small fraction c of roughness 
elements of size a with a local no-shear condition on an otherwise perfectly slipping surface. Equating the 
viscous force associated with the disturbance flow created by the defects, 0{iJ,jd'^), to the local Stokes drag on 
a defect, 0{iJ,aUe), where d ~ ajc^l"^ is the typical distance between defects and Us is the fluid velocity near 
the defects, leads to an effective slip length for the surface, A — Usj'^, given by 

A^^. (13) 
c 

When c is of order one, all these length scales (a, d. A) are of the same order and results of |134| are recovered. 
Recently more rigorous results were derived in |26|. 

When the boundary condition is locally that of no-slip, roughness shifts the position of the effective surface 
into the liquid. Calculations have been made for periodic and random surfaces jlllL I127L 11401 1164| and are 
related to earlier work on the boundary conditions for porous materials jl33l I155| and the Laplace equation 

Roughness-induced dewetting. The interaction of roughness with surface energies can lead to the spon- 
taneous dcwetting of a surface and the appearance of a super-hydrophobic state, as proposed in _36, ^37 . In 
that case, roughness could increase slip by producing regions of gas-liquid interface at the solid boundary. Let 
us consider for illustration a surface S covered with a fraction c of roughness elements of height a in a liquid 
at pressure p and let us denote by r > 1 the ratio of real to apparent surface area. In that case, the change 
in free energy, AG, to dewet the apparent area (1 — c)S between the roughness elements arises from surface 
energies and work done against the liquid 

AG - r(l - c)5(7s - 7L5) + (1 - c)S{j>a + 7), (14) 

where 7 is the surface tension of the liquid, 75 that of the solid and 7^5 the liquid-solid interfacial tension. 
Consequently, using Young's law, 7 cos 0c — is ~ iLSi we see that dewetting is energetically favorable when 
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p < —7(1 + r COS 9c) /a, which, for a given value of the pressure, will occur if the surface is hydrophobic 
(cos 6c < 0) and a is small enough. The super-hydrophobic state is therefore due to a combination of geometry 
and wetting characteristics. 

This idea is related to the so-called fakir droplets IlltJL I186| and to more general drag reduction 

mechanisms found in nature using gas bubbles |24| . Trapped bubbles in rough surfaces were studied by ^ 
in the context of contact line motion and are probably responsible for the apparent slip lengths reported in 
|182l 118:^ I184j for flow over fractal surfaces, and possibly other studies as well (see also A similar 

mechanism was quantifled experimentally using trapped bubbles in rough silicon wafers |12()| (see also the 
calculations in |18Uj ') and show promise of decreasing turbulent skin- friction drag |112| . 



5.2 Dissolved gas and bubbles 

Slip depends on dissolved gas. The amount of slip has been observed experimentally to depend on the 
type and quantity of dissolved gas. It is reported in sedimentation studies jl6j that slip was not observed in 
vacuum conditions but only when the liquid sample was in contact with air. Furthermore, the study in 
showed that tetradecane saturated with CO2 lead to results consistent with no-slip but significant slip when 
saturated with argon, whereas the opposite behavior was observed for water. Similar results were reported in 
[1611 . More generally, slip results in non- wetting systems are found to depend strongly on the environment in 
which the experiment is performed j38|. 



Flow over gas: Apparent slip. The results above, together with experiments showing dependence of slip 
on the absolute value of the pressure 1161' , and spatially- varying velocity fields TGSI , hint at the possibility of 
flow over surface-attached gas pockets or bubbles (see also the discussion in 101 ). Recent results in 38 also 
point at the possibility of flow over gas pockets associated with the contamination of hydrophobic surfaces 
by nanoparticles. We also note that the group of Steve Granick reported a contamination of their previous 
ostensibly smooth mica surfaces by platinum nanoparticles jl(J4| , possibly affecting some of their experimental 

results in una ma ma una. 

The idea of a flow over a gas layer was first mentioned in '137' and revisited in '1361 as a possible explanation 
for the attraction between hydrophobic surfaces in water: the attraction could be due to the hydrodynamic 
correlated fluctuations of the gas interfaces, analogous to the Bjerknes force between two pulsating bubbles. 
Detailed theoretical considerations have shown that it would be favorable for water between two hydrophobic 
surfaces to vaporize |1(J7| . Flow of binary mixtures have also been shown to phase separate by the sole action 
of interniolecular forces 0]. 

It is clear that flow over a layer of gas will lead to an apparent slip. Since stress must be continuous at 
a liquid-gas interface, a difference of shear viscosities will lead to a difference of strain rates. If a liquid of 
viscosity ni flows over a layer of height h with viscosity fi2, the apparent slip length for the flow above is (see, 
e.g., [1721) 

^-^-1; (15) 

Mi/M2 ~ 50 for a gas- water interface. Three differences exist however between a flow over a gas layer and 
flow over a set of bubbles: (a) The gas in bubbles recirculates, which decreases the previous estimate Eq. (|15f) 
by about a factor of four; (b) No-slip regions located between the bubbles will also significantly decrease the 
apparent slip lengths |84lll01lll23l 11241 1134| (see also 1 on the effect of non-uniform slip lengths); (c) Bubbles 
are in general not flat, which decreases the previous estimates even further (similar to the effect of roughness 
on a shear flow). 

When the gas layer is in the Knudsen regime {a ^ h <^ £m), the shear stress in the liquid, 0{iJ,j), is 
balanced by a purely thermal stress in the gas, ©(pusitth), where p is the gas density, uth the thermal velocity 
uth — 0{\/kBT /m) (m is the mass of a gas molecule) and Us is the liquid velocity at the interface. Balancing 
these two contributions leads to an apparent slip length, A — Us/'j, given by |44j 

A^^, (16) 

PUth 

which is independent of h and can be as large as microns. Note that the slip length given by Eq. (|16|l increases 
with the viscosity of the liquid. 



12 



t> 



> [» 



— o > 

>0 ^ 



40 60 80 100 120 

en 



-* E> 



— ^ 



Y(s 



Figure 2: Experimental variation of the slip length, A, with the liquid-solid contact angle, 6 (left), and the 
typical experimental shear rate, 7 (right), for the experimental results summarized in Tables [TBI Pressure- 
driven flow (o), sedimentation (•), fluorescence recovery (□), PIV (►), streaming potential (▼), fluorescence 
cross-correlations (•^), SPA (■), and APM (>). When a solid line is drawn, the experimental results are given 
for a range of contact angles and/or shear rates. Furthermore, when the value of the contact angle is unknown, 
the results are not reported. 



Nanobubbles in polar liquids? Over the last four years, many groups have reported experimental obser- 
vation of nanobubbles against hydrophobic surfaces in water 1511711176, 106. ,TT^[TMd52. 165, 166. 19Tl[Dn| . 
with typical sizes ^ 10 — 100 nm and large surface coverage (see also the reflectivity measurements in [8511143) ). 
The nanobubbles disappear when the liquid is degassed. Similar bubbles could be responsible for slip mea- 
surement in some of the experiments to date (see also jlTghl87j . How could the formation of such bubbles be 
explained? Thermal fluctuations lead to bubble sizes a ~ ^yksT/j, which are on the order of the molecular 
length. It has been proposed that shear might induce bubbles, but the mechanism is not clear ^5. An alter- 
native scenario could be a local decrease in pressure near hydrophobic surfaces due to intermolecular forces 
(see, e.g., Eq. 4.23 in gHl). 

The second important issue for nanobubbles is their stability against dissolution. A spherical gas bubble 
of radius a, diffuses into the liquid on a timescale (5 It .105. , 



DcoRT V 7 

where M is the molar mass, po the far-field pressure, D the diffusion coefficient of the gas in the liquid, T 
the temperature, cq the saturated gas concentration in the liquid (mass per unit volume) at pressure po, 7 
the liquid surface tension and R the absolute gas constant; note that this estimate is independent of surface 
tension for sufficiently small bubbles. For a 10 nm bubble, r « 10 /xs and t becomes a few hours when a « 100 
/im. It has therefore been argued that the existence of such small bubbles can only be explained if the liquid 
is supersaturated with gas In many pressure-driven flow experiments at small scale, a high-pressure gas 
in contact with the liquid is used to induce the flow; for example pressured gas at 10 atm is used to drive 
motion in equivalent to the internal pressure of a 100 nm bubble. However, the resulting equilibrium 
between a gas bubble and a supersaturated solution is well-known to be unstable as any perturbation either 
grows without limit or dissolves away. A possible resolution to the stability problem of such bubbles might 
come from intermolecular forces in the gas which become important when bubbles reach small radii and large 
pressures jl85) . 
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5.3 Wetting 



Slip depends on wetting properties. It was recognized early that friction at the liquid-sohd boundary 
should be a function of the physicochemical nature of both the solid and the liquid In particular, the 
wetting properties have been found to play a crucial role in many experiments. Wetting of solids by liquids is 
reviewed in |43l I45j and is quantified by the spreading coefficient, S' = 75 — 7 — JlSi which is the difference in 
surface energy between a dry solid surface and the same surface wet by a liquid layer (75, 7 and 7^5 are the 
solid, liquid and liquid-solid interfacial energies). When S* > 0, the solid is completely wet by the liquid and 
when S < the wetting is partial. In the latter case, a small liquid droplet on the solid surface would take 
the shape of a spherical cap with a contact line at an angle 9c to the solid, where 9c is the equilibrium contact 
angle, and is give by Young's law, 7 cos 0c = 7s — 7ls- The surface is said to be hydrophobic if 9c > 90°, and 
in that case the nucleation of small bubbles in the liquid should occur preferentially on the surface. 

Slip has been measured for systems in complete wetting |EI ^1 EZl IMl H^'*^! I126| and partial wetting 

piiniisiisiiEiEsiMiEHiEiEniinziiHziiHHiEi CHi im cnni nni nnsi ini nni iniiinsi- The amount 

of slip has usually been found to increase with contact angle, either systematically (e.g., |192| ) or only for 
non-polar liquids (5T) . All these results are summarized in the plot displayed in Fig. |2] (left), which shows 
however, overall, a poor correlation between slip and contact angle. 



The Tolstoi theory. It appears that Tolstoi was the first to try to quantify the importance of surface 
energies on slip at the molecular level ^1 153L I159| . Using concepts from macroscopic thermodynamics at 
the molecular scale, Tolstoi considered the relation between surface energies and molecular mobility (hence 
diffusivity) near a solid surface by calculating the work it takes for molecules to make room for themselves in 
the liquid, and how that changes near a boundary. The molecular diffusivity D is given by the product of the 
molecular scale, a, and a velocity, D ~ aV. The molecular velocity is V ^ a/r, where r is the typical time 
scale for hopping from one molecular position to the other. This is typically an inverse molecular frequency 
corrected for the energy it takes to create a void of size cr, which is similar to the surface energy 70-^. Near 
a solid, this energy involves the solid and liquid/solid interfacial energies, hence the possibility of having a 
different molecular mobility close to a surface. In the case of complete wetting, the Tolstoi model leads to 
the no-slip boundary condition within ± one molecular layer in the liquid, but in the case of partial wetting, 
molecules near the surfaces are found to have larger mobilities, leading to a slip length ^| I159| 

A f aa^jjl - cos9c) \ , . 

^^^"H J"' ^ ^ 

where a is a dimensionless geometrical parameter of order one, 7 is the liquid surface tension and 9c the 
equilibrium contact angle. The estimate, Eq. (|18|l . increases with the contact angle and can be orders of 
magnitude above the molecular length. 



Intrinsic (molecular) sUp. Another theory at the molecular scale uses the fluctuation-dissipation theorem 
and Green-Kubo relations to derive slip lengths from equilibrium thermodynamics [5l ll4l[T5j . Using Onsager's 
hypothesis of linear regression of fluctuations, i.e., that small fluctuations around equilibrium can be described 
by the same equations that describe the relaxation from non-equilibrium, leads to a formula for the time- 
dependent momentum correlation function in the liquid, function of both the slip length and the wall position. 
The boundary condition is found to be applied about one molecular layer inside the liquid and the slip length 
is given by 

- - (19) 

where D* = D\\/Dq, D\\ is the collective molecular diffusion coefficient, Dq the bulk diffusivity, St the structure 
factor for first molecular layer (both D* and St are dimensionless numbers of order unity), pc the fluid density 
at the first molecular layer (number per unit volume), and cls the dimensionless liquid-solid coefficient of the 
Lennard- Jones potential, Eq. Ullfl . In the case of complete wetting, the slip length is essentially zero as soon 
as the roughness is a few percent of the molecular size, but in a non-wetting situation the slip length can be 
up to two orders of magnitude above molecular size and increases with contact angle [H]- As the macroscopic 
contact angle goes to 180°, the slip length diverges as A/cr 1/{tt ~ 9cY. The theoretical predictions are 



14 



found to agree very well with MD simulations [HI C] and the approach was extended to polymer solutions in 
|128j . The details of the molecular slip mechanism were subsequently studied in llOHj which showed that slip 
at low shear rates occurs by localized defect propagation and switches to slipping of whole molecular layers 
for larger rates. 



5.4 Shear rate 

Slip depends on shear rates. In many investigations, slip was observed to depend on the shear rate 
at which the experiment or simulation was performed (see Fig. |21 right). When that is the case, the slip 
boundary condition, Eq. lO, becomes nonlinear A = A(7), and we refer to this situation as NL in Tables ITEI 
Shear-dependent slip was reported experimentally in jl2, j32, ^ £2, .116, 192. .193. .194ii I195| , with the 
strongest dependence to date in |192| . When such a dependence is not observed, Eq. ^ is linear (L) and the 
slip length is a property of the liquid-solid pair. Two drainage experiments have also reported linear boundary 
conditions, with force proportional to velocity in Eq. ((S)), but with a profile, for small separations between the 
surfaces, which differs from that given by a uniform slip length model 139). 

Most MD simulations report boundary conditions in the linear regime, except in [661 1114j . where the 
magnitude of slip was found to depend on the magnitude of the driving force, and in |158j where simulations 
give slip lengths that diverge at high shear rates A/a ^ (1 — j/'jc)'^^^- 



The leaking mattress model. A mechanical model for shear-dependent slip in drainage experiments was 
proposed in |l()()j , based on the assumption that a layer of gas bubbles is present on the solid surface. Since 
drainage experiments are unsteady, the bubble sizes will be a function of time in response to pressure variations 
in the liquid (by the combination of compression and diffusion), hereby modifying the amount of liquid which 
is necessary to drain out at each cycle of the oscillation, and therefore modifying the viscous force on the 
oscillating surface. This idea of flow over a "leaking mattress" leads to a frequency-dependent decrease in the 
viscous force, f*{w) in Eq. jSJ, given by 

/*(^) 1 _ nall{9c)coy/iJK _ cphp / ^ Coo{D - 2feo) \ 

/slip ^ , f^^ , (tj(5fc2)^y ^ npopoDa ' ^ irc^poDa \ po^o 
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where /sUp is the slip factor due to flow over bubbles, given by Eq. © {i.e., /sUp is the zero-frequency force 
decrease, associated with a slip length. A, which describes the effective resistance of the covered solid surface), 
and the other terms quantify the dynamic response of bubbles and contribute to an additional force decrease. 
In Eq. H20|l. a is the curvature of the surface, Coo is the far field dissolved gas concentration (mass per unit 
volume), po is the far-field liquid pressure, cq the dissolved gas concentration in equilibrium with gas at 
pressure po, po the gas density at pressure po, uj the frequency of oscillation of the drainage experiment, D 
the minimum distance between the two surfaces, ^ the liquid shear viscosity, n the number of bubbles per 
unit area, oq the equilibrium radius of curvature of the bubbles, k the diffusivity of the gas in the liquid, I 
a geometrical shape factor of order unity, /iq the mean bubble height on the surface and S — 1271 ^a^ fsUp / D . 
Note that f*{uj) increases with the liquid viscosity and the curvature of the surfaces. The results of this model 
compare well with the experiments of |192| . 

The critical shear-rate model. An empirical model for shear-dependence inspired by the data in |192| 
was also proposed in |146| . Slip is assumed to occur locally with a constant slip length A as soon as the local 
shear rate reaches a critical value 7c; below this critical value, the no-slip boundary condition is assumed to 
remain valid. This model has therefore slip confined to an annular region around the narrow gap where shear 
rates are the highest. With the two fitting parameters (A, 7c), the model can reconcile various experimental 
data, except for small separations of the surfaces. 



Viscous heating. In a steady flow, the rate of dissipation of mechanical energy (which, for a Newtonian 
fluid, depends on the square of the shear rate) is equal to the rate of change of internal energy due to changes 
of temperature. Since the viscosity depends on temperature, high shear rates could lead to the possibility 
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of viscous heating a flow- induced reduction in viscosity, which could be interpreted as an apparent slip. 
Assuming a traditional exponential law for the viscosity /i = exp[— /3(T — ro)/ro] and flow in a circular 
capillary of radius a, the apparent slip length due to viscous heating would be | 1U1| 




) 



{jar 



(21) 



where Tq is the reference temperature, (3 a dimensionless coefficient of order one, v the fluid kinematic viscosity, 
kt the fluid thermal diffusivity and Cp the specific heat. Although viscous heating can be neglected in most 
experiments to date, it has the potential of becoming important at higher shear rates (see also |167| on the 
issue of temperature variations). 

5.5 Electrical properties 

Apparent slip depends on ionic strength and polarity. When probing slip in electrolyte solutions 
and polar liquids, the amount of slip was found to vary with electrical properties. Sedimentation experiments 
reported that slip was only observed for polar liquids Fluorescence-correlation measurements reported 

slip lengths of the same order as the screening length k~^, which decrease with the ionic strength of the solution 
|l()8j . Drainage experiments also reported that when liquids are polar, slip fails to increase with hydrophobicity 
but increases with the dipolar moment of the liquid 31 ; this result was interpreted as a disruption by the 
drainage flow of the local liquid cohesive energy arising from dipole-image dipole interactions close to the 
surface [HJ . Finally, we note that the morphology of nanobubbles has also been observed to depend on pH 
|l(if)j (see also the discussion in 22 ). 

Electrostatic-induced averaging. When using small tracers to probe the fluid velocity, electrical effects 
need to be carefully taken into account. The flrst issue concerns measurement close to the surface. If the 
surface and the particles are similarly charged, particles will be repelled electrostatically and will not come 
within a distance ~ of the surface (if charges are opposite, tracers will stick to the the surface). This 
effect will therefore increase the average velocity of the particles when compared to what would be expected 
if electrical effects were not considered. If the averaging window has a height h > above the surface, then 
the mean flow velocity will appear increased by a factor (1 -I- I/kH), which, if interpreted as a slip length would 
give an apparent value 



Apparent slip due to charged tracer. The other potential problem with experimental methods using 
small tracers is the influence of the streaming potential on their motion 99 . If the particles are charged, their 
velocity will also include an electrophoretic component in response to the flow-induced potential difference; 
moreover, if they are sufficiently charged, this velocity will be able to overcome the electro-osmotic back-flow 
in the bulk and the particles will move faster than the local liquid, leading to an apparent slip length. If we 
consider the case of a pressure driven flow between two parallel plates separated by a distance 2/i, the resulting 
apparent slip length is given by |99j 



where Cw is the wall zeta-potential, the particle zeta-potential, q a dimensionless factor of order one de- 
pending on the ratio of the screening length to the particle size, fj, the shear viscosity of the liquid, Ue the 
electrical conductivity of the liquid and its dielectric constant. In the case of low-conductivity electrolytes, 
such apparent slip length can be as large as hundreds of nanometers. 

5.6 Pressure 

Apparent slip depends on pressure. In the velocimetry experiments of |161j . the measured slip length 
was found to decrease with the value of the absolute pressure; for water, the no-slip boundary condition 



A CwjqCp - Cti;)(er£o)^ 



(23) 
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was recovered when the absolute pressure reached 6 atm. Such results are another hint at the likely role of 
surface-attached bubbles, presumably decreasing in size with an increase in pressure. 



Slip due to pressure gradients. The possibility of surface slip due to gradients in liquid pressure was 
proposed in |137| using arguments from equilibrium thermodynamics. The idea is that the chemical potential 
of a liquid molecule depends on pressure, so a pressure gradient leads to a gradient in chemical potential, 
hence a net force F on the liquid. Assuming F ^ ksT/a, a molecular model of diffusion under force allows 
us to get the net surface velocity and estimate the slip length. For a circular pipe of radius a, the slip length 
is predicted to be 

A (24) 

where Ds is the diffusion coefhcient for molecules close to the surface and p the molecular density of the liquid 
(number per unit volume). For regular liquids such as water, the result, Eq. leads to molecular size slip 
length and suggests that the only way to get larger slip lengths would be for liquid molecules close to the 
surface to slip over a gas gap | 137| . 



6 Perspective 

Because of the great advances in micro- and nano- fabrication technologies, the ability to engineer slip could 
have dramatic influences on flow since the viscous dominated motion can lead to large pressure drops and 
large axial dispersion. As was shown in this chapter, the small-scale interactions between a liquid and a solid 
leads to extremely rich possibilities for slip behavior, with dependence on factors such as wetting conditions, 
shear rate, pressure, surface charge, surface roughness and dissolved gas. 

We conclude this chapter by presenting a perspective summarizing the interpretation and use of the slip 
boundary condition (Eq. |2Jl to describe the motion at a liquid-solid boundary. 

• Physically, there is a difference between three different types of slip: (a) Microscopic slip at the scale 
of individual molecules, (b) actual continuum slip at a liquid-solid boundary {i.e., beyond a few molec- 
ular layers) and (c) apparent (and effective) slip due to the motion over complex and heterogeneous 
boundaries. 

• From a practical standpoint however, the distinction is not important. Whether it is real slip or apparent 
slip due to the interplay of many physical parameters, we have seen in this chapter that a large number 
of (generally small) experimental systems display some form of reduced resistance to fluid motion. 

• The (apparent) slip lengths reported experimentally span many orders of magnitude, from molecular 
lengths up to hundreds of nanometers. The impact of slip on systems with typical dimensions larger 
than tens of microns will therefore likely be limited. 

• Molecular theories are able to predict intrinsic slip lengths of up to tens of nanometers for hydrophobic 
systems, suggesting that any measurement of larger slip is affected by factors other than purely fluid 
dynamical. 

• The parameters that contribute to apparent slip include roughness-induced dewetting, the amount and 
nature of dissolved gas, contamination by impurities and viscous heating. Other parameters that influ- 
ence the magnitude of apparent slip include contact angle, shear rate, electrical properties and pressure. 

• Finally, although it is usually assumed that slip only occurs on hydrophobic surfaces, a large variety of 
hydrophilic surfaces with different wetting properties have be shown to be prone to slip (Fig.|21). 

Other more complex behaviors remain to be understood, including dependence of the results on the molec- 
ular shape and size [501 1551 1195j . probe size |108| . or viscosity (lUl I116| . The development of alternative 
direct experimental methods would allow for a more precise quantification of slip phenomena. Similarly, it 
might be valuable to reproduce some of the experiments discussed above in degassed and clean environments 
to quantify the influence of dissolved gas on apparent slip. Answers to these questions will probably allow for 
a precise engineering of slip in small-scale systems. 
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